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1 Introduction 



Two dimensional quantum electrodynamics or the Schwinger model may be served as a 
laboratory for studying four dimensional gauge theories and important phenomena such as 
confinement, screening, chiral symmetries and • ■ ■. On a flat space-time, it is known that in 
the massless Schwinger model via a peculiar two dimensional Higgs phenomenon, the photon 
becomes massive and the Coulomb force is replaced by a finite range force, giving rise to 
the screening phase. When the dynamical fermions are massive, the model is not exactly 
soluble, but a semiclassical analysis reveals a linear rise energy for opposites test charges q 
and q, binding them into qq pairs 0. By taking into accounts the finite distance corrections, 
one can show that besides the linear confining term, the potential is also composed of a 
screening term which is modified with respect to the massless case ||. 

One of the interesting questions is how a curved background can modify these effects? 
We think that this is an important question, because it can be viewed as a first step in 
studying these physical effects in the context of quantum gravity. Moreover, they may have 
applications in string theory and quantum gravity coupled to nonconformal matter (note 
that the kinetic term of the gauge fields spoils the conformal invariance of the theory). 

Although the Schwinger model on curved space-time has been studied in several papers 
but the confining aspects of the model is still unclear (for a discussion about the subtleties in 
determining the confining phase of the Schwinger model on curved space see For example 
the author of || has suggested that the curvature of the space does not change the confining 
behavior of the massive Schwinger model, and in || it has been claimed that on a particular 
black hole, the massless Schwinger model remains in screening phase. Also in |4j], by com- 
parison the role of temperature and the curvature, it has been argued that the curvature 
may modify the confining or screening nature of the Schwinger model. More recently, this 
model has been investigated on a constant negative curvature space-time (i.e., the Poincare 
half-plane) 0, and it has been found that the confining feature of the Schwinger model 
depends on the geometry of the space-time, for example, in different regions of Poincare 
half-plane the system is in different phases (i.e. confining or screening phases), which is in 
contradiction with flat space-time results. Also in |§, the confining aspects of the mass- 
less Schwinger model on de Sitter space in different coordinates has been discussed and the 
results have been extended to non-abelian situation. 

In this paper we want to study the massive Schwinger model on a general noncompact 
Riemann surface, and try to obtain, as much as possible, information about the confinement 
in this background. We consider only the effects of geometry on confinement phenomenon 
and ignore nontrivial topologies. 

In the Schwinger model the potential of external charges can be obtained using several 
different methods, for example : (i) Integrating over the fermionic fields (or equivalently 
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over bosonic fields in the bosonized version of the model), results an effective action for 
the gauge fields, from which the potential can be extracted by solving the corresponding 
equations of motion, (ii) Integrating over gauge fields in the bosonized action and then 
obtaining the static solutions of the equation of motion and deducing the potential of the 
system as the difference between the Hamiltonian in the presence and absence of external 
charges respectively. In this paper we follow the second procedure. 

The plan of the paper is as follows: In section 2, our main task is to obtain the bosonic 
version of the massive Schwinger model on a general non-compact curved background. To 
find this bosonic representation, it is necessary to employ an appropriate normal ordering 
description in defining the composite fields, which are always present in the bosonic repre- 
sentations. As we will see, this needs several change of field variables in order to reduce our 
main model (interacting massive fermionic fields in curved background) to a free massless 
fermionic field theory in flat space-time. Then using the equivalence of this theory with a 
massless bosonic field theory, we may obtain the appropriate normal ordering description, 
from which our bosonic representation can be found. 

In section 3, by introducing two external charges q and q in the bosonized action, and 
by restricting ourselves to static metrics (so that the static potential between probe charges 
becomes meaningful), we obtain the energy of widely separated external charges for small 
dynamical fermions mass. Obviously, the final result depends on the metric. For a specific 
example (the full de Sitter space-time), we will show that one recovers a confining behavior 
in the above mentioned limit. We also show that when e'/q G Z, where e' (q) is the charge 
of external (dynamical) charges, the external charges can not modify the energy .This is the 
same result as the flat case. Note that in section 2 where we use the path integral approach, 
it is more convenient to consider the Euclidean signature, but in section 3, where we are 
going to calculate the energy, it is better to work with Minkowskian signature, in which the 
concept of energy is more natural. 

To study the same problem but with external charges with finite distance separation, in 
section 4 we restrict ourselves to the case in which the gauge coupling is very large with 
respect to the variation of the metric. The main reason for such a choice is that, the massive 
scalar Green function, which is needed in calculating the energy, is only known for a few 
number of space-times. But if we restrict ourselves to the above mentioned metrics, we can 
use the WKB approximation to find a general expression for the energy of external charges. 
Like the flat case, the energy takes a screening correction terms beside the confining term. 
As a result we show that at m = 0, the phase structure of the system depends on the metric 
and in contrast to the flat case, the model can be in confining phase in specific cases (such 
as the full de Sitter space-time) which is in agreement with the result of ||. 
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2 Bosonization of the massive Schwinger model 



All two dimensional spaces are conformally flat, hence any noncompact Riemann surface is 
described by the metric 

ds 2 = ^g(dt 2 + dx 2 ), (1) 

where y/g is the conformal factor. On this space-time, the massive Schwinger model is 
described by the partition function 

Z = J DApDipDipi exp{— J d 2 x[^{ftY{d^ - iqAJip + mftft) + ^^Ih ( 2 ) 

where 7 M = 7°e^ are the curved space counterparts of Hermitian Dirac gamma matrices: 
7° = cr 2 ; J 1 — (J\. <Ji are Pauli matrices. <9 M is the covariant derivative and the zweibeins 
are defined through 

9r = effiat, ST = <^ ah - (3) 

The metric components are g^ v = (1 / ■ s /g)S tlu , where 5 tJ,u = 5^ u = diag(l, 1). q and m are the 
charge and the mass of dynamical fermions respectively. In eq.(^), the dimension of q and 
m is the inverse of the length dimension and e is dimensionless. The (dual) field strength 
F, is described through F = e^d^Ay in terms of gauge fields A^, where e^ u = e M „ and 
e 01 = — eio = 1. The partition function (0) can be written as 



°° (—m\ k k r i 

Z = Y. [ -JT < II J fxj^MtfteMxj) >, (4) 



(—to)* 
k=o h - j=l 

where x = (x,t). In eq.(^), the expectation values are computed from the Lagrangian 



L = (fy - iqA^ + ^=_F 2 . (5) 

In terms of the new fermionic variables ip an d ^ ■ 

i> = exp(g7 5 (^)^, 

ft = ft exp(q>y 5 (p), (6) 
where A^ = e^ v d u ip and 75 = — ry 1 ^ , the Lagrangian @ can be rewritten as || 

L = Jgft-fd^ + ^g(-^ipA(p + ~<pAA<p), (7) 
in which A is the Laplace-Beltrami operator: A = g^ u d^d v and 



3 



In this way the massless part has been became a free field theory on a curved background 
with an effective action containing an anomalous term (— (^l 2 /2)^gipA(p) coming from the 
Jacobian of the transformation . The term (l/2)y^^A(A — fi 2 )<p, is the effective La- 
grangian density of the gauge fields. The fermionic part of the Lagrangian fl7|) is free and 
hence is invariant under Weyl transformation: g^ u — > VPg^ and ip — > Vt~^ip. Choosing 
Vt = g~4, one obtains 

4> — > A = g^ijj, ^ — > A 1 " = g*ffi. (9) 
Therefore the partition function (fj) can be reduced to 



E -^T~ <UJ d 2 x^{^)\\^) exp[-2g 75 ^(x i )]A(x i ) >, (10) 



(— m) 
fc=o i=i 

where the expectation values are computed from 



L = \^ a d a \ + i^A(A - /i 2 V. (11) 

Therefore the massless fermionic part of the theory is reduced to a free fermionic field theory 
on a flat space-time. Using (7s) 2n = 1, it can be easily shown that 

A t exp(-2g 75 ^)A = A+[^-^ exp(-2^) + exp(2g^)]A. (12) 

The fermionic part of ( [□] ) is chirally invariant, hence the only non-zero terms in ( jTO]) are 
those with equal number of a + and ct_ |jj~0f , where a± = A^(l ± 7 5 )A/2. Hence 

z = E Tina < II / rf S rf2 %^^( x i)^^(yi) ex p[ 2 ^(yi)] ex Ph 2 ^( x j)] CT +( x j) CT -(yj) > • 

fc=0 v "/ j=l 

(13) 

The expectation values of the fermionic part can be established in the same manner as in 
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/ D\D\1 n Mx>-(y.) ex P (- / d 2 xX^d a X) = (1) ^^' ^ y > ? . (14) 

J j=i J Z7r iiij=ii x i yjj 

But the Green function of the massless scalar field </>(x) is < 0(x)0(y) >= —1/ (47r) ln[e 2 (x — 
y) 2 ], where the mass e 2 (e 2 — > 0) is introduced to avoid infrared divergences. This is 
equivalent to infrared renormalization of exponential of massless scalar field discussed in [|12| 
in operator language. So using 

/a 2 I r 

£>0(x) exp[£j =1 -^£>(xj, x,)] exp[Ej =1 ia^( Xi )] exp[-- y (9^(x)) 2 rf 2 x] = 
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^ j=1 3 i<j 



(15) 



in which we have used e 2 — > limit and -D(x,, Xj) = — 1/ (47r) ln(xj — x^) 2 , one can write the 
eq.(|l4|) in the form 

(_L)2fc eX p[27rSf =1J D(x l , x*)] exp[27rSf =1J D( yi , y,)] x 



D0(x)exp{2z v ^E(0(x J )-0(y J ))]}exp[- / d 2 x-(^(x)) 2 
j= i J 1 



(16) 



Hence the eg. (|i~3|) can be written as 



m 



2k 



< n / ^ rf2 %^*( x j)^*(yj) ex p[ 2 ^^( x i> x i)] ex p[ 2ff ^(yi>yi)] 



S(27r)»(*!) J=r 
exp[-2g^(x i )] exp[2g<p(y_j)] exp[2z v /7T0(x i )] exp[-2iv / 7r0(y J )] >, 



x 



(17) 



where the expectation values is calculated from the following Lagrangian (note that ^g^ u = 

L = ^^A(A - /i> + hjggTd^W. (18) 
Now if we note that there is the charge conservation law in contracting the composite 



fields e 1 ^ ||10| , then the partition function (|17D can be written as 



rn 



2k 



£ (27T)^(fc!)2 



< if / ^•'•// 2 //,/'(X ; )//(>'; 



X 



exp[2iVvr(0 + — y>)(x i )] exp[-2iV7r(0 + — <p)(yj)] > 

JlT JlT 



< E tt{- / ^( x ) cos[2v^(0(x) + -^(x)]} fc 

< exp{— / d 2 xh(x) cos[2v/7r(0(x) + -j=(p(x))]} >, 

IT J \ TT 



> 



in which 



(19) 



(20) 



/i(x) = <7 4 (x) exp[27T-D(x, x)]. 

By changing the field variable — > — (iq/\/W)<p, and using = e^ u d u ip, the partition 
function ([TOD becomes (note that this change of field variable is somehow inverse of @ and 
couples the matter and gauge fields) 



Z = J DA^Dtexpi- J d 2 x^-g^d^dA - i-^e^ A„d v (f> + ^-F 



vn i 

-— g _3 (x) exp[27rD(x, x)]cos(2 v / 7T0(x)]}, 

TT 



(21) 
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where & v = {\j ■ s fg)e ia ' . By integrating over the gauge fields one can see that the effective 
Lagrangian of the field is 

Left. = \^gg^d^d v (t> + ~ ~9* exp[27r£>(x, x)] cos[2^(x)]. (22) 

Now we want to convert the mass term in eq.(^) to a more convenient form, by absorbing 
covariantly the ultraviolet divergence in the mass m or in operator language, by defining an 
appropriate normal ordering prescription (N^) with respect to the scale /J |T3J. To do so in 
the flat case, we write exp[27i\D(x, x)] as 

exp[27rD(x, x)] = exp[— 27r(G M (x, x) — -D(x, x))] exp[27rG A1 (x, x)] 

= exp( 7 + ln|)exp[27r^(x,x)], (23) 

where is the Green function of a massive scalar field with mass /x; which for ^fg = 1 is 
^fflat _ i / (27r)i^ (A''l x _ y|)- Kq is the modified Bessel function of the second kind and 7 is 
the Euler constant. We can now absorb exp[27rG M (x, x)] in redefinition (renormalization) of 
mass m: 

mV>V = M g6X ^ 7 cos(2v^F0), (24) 

27T2 

where M = mexp[27rG^ at (x, x)]. In operator language |IJ, this is equivalent to the well 
known result j0|: 

m^V(x) = ~m g 6XP f 7 ^ Nfj, cos[2 Vtt0(x)] . (25) 

27T2 

Now let us consider the normal ordering in curved background. In this case, one can use 
a covariant point-splitting method in order to regularize the composite operator <p 2 (our 
regularization method should respect general covariance) 

< r 2 eg .(x) >= Hm [G„(x, x') - G™ (x, x')], (26) 

where G^ s (x., x') (extracted from DeWitt-Schwinger expansion) is the counterterm needed 
to regularize the ultraviolet divergence of the Green function [|13|] 

G™ (x, x') = -i-[2 7 + ln(/A 2 ) - + 0{e% (27) 

R is the scalar curvature of the space and e is one half of the proper distance between x and 
x'. ( Note that in the flat case, R = 0, eq.(p7|) coincides with linv^x G^ at (x, x'), and from 
(|26| ) we obtain the usual regularization which kills out the loops.) Now if we follow the same 
steps as in the flat case, but here using G^ s (x, x) instead of G^ at (x, x) , we obtain 

m^(x) = -m—^— exp{-27r[Gj 5 (x,x) - D(x, x)]}^ cos[2V^>(x)], (28) 

7T<7*(X) 
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in which we have defined the normal ordering as 

iV M cos[2v / ^0(x)] = exp[27rGj )5 (x,x)] cos [2^0 (x)]. (29) 

In the static curved space-time we can use a normal ordering which like the flat case kills 
the loops (the vacuum is now defined using the global time-like killing vector of the static 
space—time ||14|| ). In the same manner as the flat case, we obtain 

mipty = -exp[27i(D - G M )(x, x)]^ cos[2^(x)]. (30) 



It is easy to show that (|30|) is equal to (g|). From (^) (or equivalently (|28| ) for static 
space-time), we obtain 

< ^V(x) > m=0 = exp[-27rG(x, x)], (31) 

where G(x, x) = G /1 (x, x) — D(x, x). So, as we expect, the expectation value < if)ty > is 
independent of the method of regularization ( or normalization) of composite fields. 

3 Quark— antiquark potential in the massive Schwinger 
model on static curved spaces: widely separated charges 

In this section we obtain the energy of two external static charges introduced into the massive 
Schwinger model. We consider an infinite static conformally flat space-time, with trivial 
topology, described by the metric 



ds 2 = yfg(x)(dt 2 - dx 2 ). (32) 

The bosonized action of the massive Schwinger model, in the presence of the covariantly 
conserved external current 

e' 

J°(x) = — (-5(x-b) + 6(x-a)), J 1 = 0, (33) 
V9 

describing two opposite point charges — e' and e' located at x = b and x = a respectively, is 
d 2 x[-Jg-g^d^d u <P+ ^=F0+ -g% exp[-2vrG(x, a;)]A^ cos(20F<Kx)) 



2^ 

where 



^-^—F 2 + fjF}, (34) 



ij(x) = ef[e(x-a)-e(x-b)) = {^ ( 35 ) 

[0, xf[a,b\. 



Eq.(^4j) is the Minkowskian version of ( pi]) in the presence of external current (|33|). In 
deriving ([34]) we have also used eq. (|30|) . Integrating over the field F results 



/ ] <Yj\ i 

d 2 i[-Vjf^9^ + exp[-2vrG'(x, x)]^ cos(2v^0(x)) 
2 7T 



2tt 



fa + 01, (36) 



in which we have defined rj = (y/ir/q)?}. For widely separated charges, fj is equal to the 
constant e' in the whole space (in the next section we will obtain the necessary corrections 
for finite charge separation distance). If one changes the variable — > <fi — rj, finds 

S'eff. = / rf 2 x[^^0^0+-( 7 3exp[-27rG(x,x)]iV / ,cos[2 v ^(0-r ? )(x)]- 

J Z 71 

yv^ 2 ]- (37) 

For (e'/q) G Z, the action is not modified by the presence of external charges, hence the 
energy is not changed. In other words, in this case, external probes with charge e! are 
screened by dynamical fermions. The same effect occurs in flat space-time @. 

Now let us restrict ourselves to |0| <C 1 scalar fields. In this regime, the action (|36|) 
becomes Gaussian and the classical solutions of the action coincide with the quantum ones. 
The classical equation of motion for static field <p is 

[d\ — n 2 ^fg — 4v / vt7bE- v /(? cos(2\/vr?7)]0 = — 2v / 7rm£ v /(7 sin(2\/7n7), (38) 

in which E is defined through 

E(x) = -j- exp[-27rG(x, x)\. (39) 
g*7i 

The solution of eq.(^) is 

= 0" 1 2 v / ^mE v ^sin(2 v / ^r/), (40) 

where 

O = -df + {i 2 y/g + Amn^^g cos^^r/). (41) 
Now as the eq . ([40]) can be written as 

(f) = [1 - 0~ l {0 - 2^mZ^)} sin(2v/^), (42) 

therefore the condition |0| <C 1 is satisfied when mE <C fi 2 . The energy of the system is 

E = [ T°dx = - f Ldx, (43) 
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where T£ is the energy momentum tensor and E is the energy measured by a static observer 
with respect to the coordinate fl32[) . By substituting fl40|) into (|36|) , we obtain 

rb _ 

-£?(//) — — I [2nm 2 sin 2 (2^Fr/)S v ^O" 1 S v ^ + mE v ^cos(2v / 7r^)]rfa;. (44) 

J a 

Up to the first order of m and by considering mE -C /x 2 , one can arrive at 

rb 

E{rj) = -mcos(2 x /7rr/) / E(x)v^dx. (45) 

J a 

Therefore the energy of external charges is 

E ext . = E(rj)-E(0) 

= m[l - cos(2tt-)] f E{x)y/gdx. (46) 

Q J a 

This completes our general result for the energy of widely separated external charges ±e' 
in a static curved background. This result coincides exactly with the result of [Q, in which 
the same problem has been discussed using the fermionic action. As is clear from (fl^), the 
determination of the energy is now complicated by the presence of the conformal factor (y/g) 
and S, which the latter can be expressed in terms of Seeley DeWitt coefficients although not 
explicitly known for general curved space-times. To obtain an insight about this result, let 
us consider a specific example. 

Example: complete de Sitter space— time 



Consider the following geodesically complete space-time [[16 



which represents the full two-dimensional de Sitter space-time, with scalar curvature R = 
2/A 2 . By the change of coordinate x — > x(r), where dr 2 = dx 2 / (1 + x 2 / X 2 ) 2 or r = 
A cos~ 1 [A/(a; 2 + A 2 ) 2 ], this space takes the conformally flat form 

2 

ds 2 = (l + ^-)(dt 2 -dr 2 ) 

(dt 2 -dr 2 ), < — < — . (48) 



cos 2 (^r~" 2 ' A ' 2 



On spaces with constant positive curvature R, G M (x, y) is [T7 

1 D,2 -r -r 

hm G,(x, y) = -—{27 + ln(— ) + *(- + a) + *(- - a)}, (49) 
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where a 2 = 1/4 + 2fi 2 /R and \I/ is the digamma function. As a consequence £ is 

£ = 1 exp[ 7 + Iln(|) + ^(1 + a ) + U>{ 1 - - a)]. (50) 

The energy of widely separated external charges is then 

e' 

E cxt , = m[l - cos(27r-)]S(6 - a). (51) 

b and a are expressed in terms of x coordinate. Although the energy is not linear in terms 
of charge separation distance d = j^dx/ \Jl + x 2 / A 2 , but for d — ► oo we have — ► oo, hence 
the system is in confining phase. 



4 Energy of finite separated charges 



If we consider (ftq) for m = 0, it gives -E ex t. = 0, which is not right as it does not contain 
the finite separation contributions. In this section we want to study the case in which the 
external charges have finite separation, although the energy (f45D , is yet the dominant long 
range term. In this case, the function fj can not be considered as a constant in the whole of 
the space. As has been mentioned in section one, in flat space-time the energy expression, 
beside the confining term, consists also of correction terms which are due to the screening 
behavior of the system. Here we want to calculate these correction terms for a static curved 
space-time. 

Let us consider the general expression (0), in which we must now treat fj as a x-dependent 
expression. Before any calculation, we must say something about G(x, x) in (|36|). As was 
pointed earlier, we don't know the explicit form of this function for a general curved space- 
time. So we will restrict ourselves to a specific condition for the metric (or the gauge 
coupling). We first note that on the conformally flat space-time ([[]), ^^(x, x') satisfies 

(-^2 + ^2- /^)G,(x,x') = -5(x,x'). (52) 

Using the expansion 

G M (x, x') = i- J G^k, x, x')e ik ^dk, (53) 
we obtain the following equation for G M (fc, x, x') 

[f_ _ ^ + ^^g)]G,(k, x.x>) = -5(x, x'). (54) 
Now if we consider the large gauge coupling regime, such that 

dq^ , \ 

dx~ ^ ^ (55) 
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then the solution of ([54]) ( at zeroth order of WKB approximation) is 

G^k,x,x') = —— —j= ! —==-exp[-\f ]Jk 2 + ^^{^du\}. (56) 

2{k 2 + fi 2 ^g(x))*(k 2 + fi 2 ^g(x'))4 Jx v 

In the limit x — > x', G^{x, x') is found to be 

hm G,(x,x') = i- lim K {^[(x - x'f + (t - t'f]}^. (57) 

x^x' Z7T (x,t)-*(x' ,t') 

where 5s = g^(x)[(x — x') 2 + (t — t') 2 ]^ is the distance between two points x and x' when 
x' — > x. Now using -D(x, x') = — (l/27r)ln(5s) and by considering the behavior of Kq for 
small arguments, G(x, x) (introduced after eq. (|31~D ) is found to be 



27rG(x,x) = -{ 7 + ln^}. (58) 

Thus eq.([Jl]) yields 

< *V (59 ) 

27T2 

Note that in this approximation, the metric factor in eq.([3lD is canceled out by the corre- 
sponding term in G(x, x), and the final result is the same as the flat case as one expects 
from the results of f\T 



Now let us go back to eq. (|36|) . If we restrict ourselves to the regime (|55|), in which eq.(p8|) 
has been obtained, and change the variable <fi as — > <fi — rj, the effective Lagrangian fl36|), 
for static fields with \<f> — r]\ <C 1 condition, is reduced to the following Lagrangian 

+ -rid\ri — 2nmT l ^/gr] 2 . (60) 
Using (|39|) and (|58"1), £ is obtained as £ = qexp("f)/(2Tc 3 ^ 2 ). The field equation is 

Using the following identity which holds at large coupling limit 



[Jg{x)} a r I — 7= 7= 

I X r - r— r^T^ exp [~ / V^ygiu) + AnmEyf g(u)du]dx 

[/i 2 Jg(x) + AnmT,^ g(x)}! 3 Ja v 



g(x)] c 



[n 2 Jg(x) + 4irm'EJg(x 



— — r ex P[- / V V?\Jg{x) + 4:irmEy/g(x) dx] + c,(62) 
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in which c is an arbitrary constant, we obtain 

^£[F(x,b)-F(x,a)], x>b, 

(j)(x) = { v^F^ - ^^[2F{x,x) -F(x,a) - F(x,b)}, a < x < b, (63) 

-jE-l£[F{x,b)-F{x,a)], x <a, 



where 



and 



/■'(,-.,) .-= v^y) H f V f\(u)du\), (64) 

f*(x)f*{y) J* 



f(x) = (fi 2 + Airmi:)^g(x). (65) 

In deriving fl63l), the Green function ( |56| ) has been used. For a fix -F(a;, y) is an increasing 
function of y in ?/ < x region and decreasing in y > x. So F(x, y) < F(x, x). But /i 2 F(x, x) = 
fi 2 /([i 2 + 47rmS) < 1, which is not necessarily small. Therefore the condition \<p — rj\ <C 1 
is satisfied only when the factor e'/q is very small. The same arguments holds when e'/q is 
close to an integer number. 

Putting (|63"D back into the Lagrangian flBTf ) gives 

E ext . = ^{[ [V9V 2 ~ J, 4 (66) 



Using (|62|) we arrive at 



,e' r/1 /i 2 ./i 2 



H — -l—o — 5 3-expf— / t*{u)du)\\. (67) 

4 fHa) /§(&) (f(a)f(b))i Ja 

The energy of external charges is composed of two parts: The first part is proportional 



to f a HyJ g(x) dx and when mE <C ft 2 , coincides with the confining term ( [45] ) (in the limit 
e' < g). On the flat space-time the remaining terms are due to the screening of external 
charges by dynamical fermions. On the curved space, the problem is more complicated and 
we can have a confining situation even in the massless Schwinger model. To see this, assume 
that m — 0. The energy of external charges is then 

Eext. = ^-[gH a )+9Hb) -2g*(a)g*(b)exp(-f* ^ g*(u)du)], (68) 

which for largely separated external charges, in contrast to the flat case, does not tend to a 
constant, and depends on the value of the metric at the position of external charges. The 
condition of validity of the eq.(|68|) is only fl5"5|). In other words, in m = there is no need to 
condition (e'/q) <C 1 (or \<f> — rj\ <C 1) in deriving fl6"8|) . For example in the massless Schwinger 
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model on de Sitter space ( fRf) , the energy of largely separated charges (in the coordinate x) 
is 

e /2 / I2" / 52" 

£ - = i^f + If + V 1 + # (69) 
Although this energy is not linear in terms of charge separation distance, but for d — > 00, E 
is infinite and the system is in confining phase. 

From eq.(^), we can conclude that in contrast to the flat case, in which the energy 
of external charges is only a linear function of charge separation distance, on the curved 
space-time the energy depends also on the position of external charges. As it is clear, the 
infinity of the energy in a confining situation is related to the separation of charges, but 
the increasing rate of the energy with distance is not unique in all regions. For example in 
the region x ~ x sing , where x sing , G R\J±oo is a point at which the metric is singular, two 
external charges located close together can have a finite energy but by moving one of the 
charges, the energy increases very rapidly. 
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